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Abstract 

The one loop quantum corrections to the classical reflection factor of the sinh- 
Gordon model are calculated partially for general boundary conditions. The 
model is studied under boundary conditions which are compatible with integra- 
bility, and in the framework of the conventional perturbation theory generalized 
to the affine Toda field theory. It is found that the general form of the related 
quantum corrections are hypergeometric functions. 
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1 Introduction 



Affine Toda field theory is an integrable quantum field theory in two-dimensional 
Minkowski space-time which possesses remarkable properties and rich algebraic struc- 
ture (for a review see |J). This theory becomes more interesting when it is 
restricted to a half-line. However, for most of the Toda theories corresponding to 
affine simply-laced algebras, the boundary conditions are limited to a finite number 
which preserve integr ability Corrigan et.al have classified the boundary condi- 
tions which preserve classical integrability. However, there still remains much to be 
studied in relation to quantum integrability on the half-line. For models based on an 
much is now known [^,|T^-|T^[. The simplest affine Toda theory, the sinh-Gordon model 
has been studied much more than other models in the context of integrable boundaries. 
This model is the only theory in the ade series of affine Toda field theory for which 
continuous boundary parameters are possible. 

In recent years there has been considerable interest [ 12 . 14 - 18] in perturbative affine 
Toda field theory. The motivation behind this fact is that the boundary S-matrices 
of the models are largely unknown. The most progress has been made for affine 
Toda field theory for which the general form of the boundary S-matrix has been found 
by Ghoshal [|19[]. In fact, the boundary bootstrap equations yield the boundary S- 



matrices up to some unknown parameters. The perturbation method not only provides 
an additional check of the results which come from the bootstrap technique, but also it 
could make a connection between the unknown parameters of the boundary S-matrices 
and the boundary parameters which are involved in the Lagrangian formulation of the 
theories. 

Firstly, Ghoshal and Zamolodchikov [|2D| obtained the soliton reflection factors in 
the sine-Gordon model with a boundary consistent with integrability. Then, Ghoshal 
[W\ using these results calculated the reflection factors of the soliton- anti-soliton bound 
states (the breathers) of the model. One of the interesting problem in the boundary 
sine(sinh)-Gordon model is to find the relation between the free parameters appearing 
in Ghoshal's formula and the boundary data appearing in the Lagrangian of the model. 



Corrigan |14[ was the first to notice that the lightest breather reflection factor of the 
sine-Gordon model is identical to the reflection factor of the sinh-Gordon model after 
an analytic continuation in the coupling constant. In a recent paper Corrigan and 



Delius [|13|] studied the boundary breather states of the sinh-Gordon model on a half- 
line. They calculated the energy spectrum of the boundary states in two ways, by 
using the bootstrap equations then by using a WKB approximation. By comparing the 
results obtained by the two methods, they provided strong evidence for a conjectured 
relationship between the boundary parameters, the bulk coupling constant and the 
parameters appearing in the quantum reflection factor calculated by Ghoshal. They 
carried out the calculations in the special case when the boundary parameters are equal 
and the boundary condition preserve the <p —<P symmetry of the bulk theory. 



In [14] the quantum corrections up to 0(f3 ) to the classical reflection factor of 



the sinh-Gordon model were found when the boundary parameters are equal. In this 
case, the static background configuration is = 0. If the boundary data are different 
then, the lowest energy solution will not be a trivial background. The correspond- 
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ing perturbation theory involves complicated coupling constants and two-point Green 
function as well. Recently |21] the quantum reflection factor has been calculated in 
one loop order up to the first order in the difference of the two boundary parameters. 
The result |2T] provide a further verification of Ghoshal's formula. This paper extends 
the results of by calculating the quantum reflection factor for any value of the 
boundary parameters. It is found that most part of the related quantum corrections to 
the classical reflection factor may be expressed in terms of hypergeometric functions. 
This result and how it could relate to Ghoshal's formula is discussed in the conclusions. 



2 Boundary sinh-Gordon model 

The sinh-Gordon theory on the half-line is a massive scalar quantum field theory in 
1+1 dimension whose corresponding untwisted affine Kac-Moody algebra is . The 
Lagrangian density of the theory is: 

£ = 6{-x)£ - 5(x)B (2.1) 

Here, £ is the bulk Lagrangian density of the model which is given by 



£ = -d^(j)-V{(j)) 
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1 2m 2 
= -d,cf>d^-— cosh(/?a0), (2.2) 

where m and f3 are a mass scale and a coupling constant of the theory. Moreover, the 
boundary potential B has the generic form |2^ 



B = J2 (*oe-> + axefr) . (2.3) 

In the above relation, the two real coefficients a and o\ are arbitrary and indicate ||[7j 
the degrees of freedom allowed at the boundary. In fact, Bowcock et.al || obtained 
some results about the form of the the boundary term via a generalized Lax pair when 
there is a boundary. For further discussion on the boundary parameters see P,[7|. 

The sinh-Gordon model is integrable classically which means there are infinitely 
many independent conserved quantities Q± s where s is an arbitrary odd integer. On 
the other hand, the model is integrable after quantizing which implies the S-matrix 
describing the n-particles scattering factorises into a product of two-particles scattering 
amplitudes. The S-matrix describing the elastic scattering of two sinh-Gordon particles 
of relative rapidity 9 is conjectured to have the form |l, 22,23 



s{0) = -mhf) (2 ' 4) 



where we use the hyperbolic building blocks 



ITTX 



(x) = Sinh(g/2 + 4 j (2 5) 

v ; sinh(0/2 - ) ' v ' ; 
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9 /3 2 

and the quantity B is related to the coupling constant (3 by B = ■ 

In order to maintain the integrability on the half-line, the boundary potential must 
satisfy the following equation 

d<f) dB , N 

or 

^ = _^i( ffie «^_ ffoe -*/^ atI = , (2 . 7) 

where we use the normalization condition a 2 = 2 which is customary in affine Toda 
field theory. In what follows the dimensional mass parameter m will be taken to 
unity. It is also convenient to use cr, = cos a^. For the boundary sinh-Gordon model, 
besides to the two-particle S-matrix it is necessary to know the boundary S-matrix or 
reflection factor describing one particle reflection off the boundary. Firstly, Ghoshal 
and Zamolodchikov [[HJ calculated the soliton reflection factors for the sine-Gordon 
model by solving the boundary Yang-Baxter equation. Then, Ghoshal [TI| calculated 
the soliton-antisoliton bound state reflection factor. He used the boundary bootstrap 



equations along with the result of reference [20|. The general form of the quantum 



reflection factor in sinh-Gordon model may be derived by regarding the lightest breather 



reflection factor of the sine-Gordon model ||19|| , calculated by Ghoshal, and performing 
ilytic c 

K q {9) 



analytic continuation in the coupling constant to find 

(l)(2-B/2)(l + S/2) 



(1 - E(a , a u (3)){l + E(a , a u (3)){l - F(a , <r x , (3)){l + F(a , a u 0)) ' 

(2.8) 

Note the bulk reflection symmetry leads to F = or E = when o"o = <J\ (note only 
one of them vanishes). In fact, the exact form of the E and F is an open and hard 
problem. Recently Corrigan and Delius |13] obtained the function E in the special case 
when cr = (J\ = cos air and \ < a < 1 as 

E = 2a(l-B/2). (2.9) 

They found the above formula by equating the results of the WKB approximation 
method and the bootstrap technique. 



3 Low order perturbation theory 

For affine Toda field theory the perturbative calculation is performed around the static 
background field configuration, so standard Feynman Rules may be used. By expanding 
the bulk and boundary potentials in terms of the coupling constant /3, the three and 
four point couplings can be deduced. We find for the sinh-Gordon theory 

C£L = ^/Mrm(VWo) (3.1) 
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^k = ^ 2 cosh(v/2/30o), (3.2) 

where <po represent the background solution to the equation of motion of the model 
and similarly 

CgL*n, = (V ^ ~ W^ 1 *) (3-3) 

Cit ndary = | + a e-^) . (3.4) 

On the other hand, the static background field can be found through linear perturbation 
of the equation of motion and the boundary condition of the model @,|14[ to obtain 

, r- 1 _|_ p 2(a;-xo) 

e Wo/V2 = £ + ^ (3 5) 

where the parameter xq is related to the boundary parameters by 



cothxo = \ —— — -• (3.6) 
V 1 + CTi 

So, the three and four point couplings corresponding to the bulk potential take the 
forms _ 

C ( £lk = -^cosh2(x - xo) (coth 2 2(x - x ) - l) , (3.7) 

^L = ^ 2 (2coth 2 2(x-x )-l). (3.8) 
In the same manner the three and four point couplings of the boundary are given by 

\/~2{3 

C boundary = (^i coth x - ct tanhx ) , (3.9) 

4 ft 2 

C ioundary = (^1 C ° th X + ^tanl^) . (3.10) 

The next step is to find the propagator for the theory. It has been shown [|14]] that 
the two-point Green function for the sinh-Gordon model on a half-line is 

+ KJ(-k,x)f(-k,x')e- ik ^ +x '^ , (3.11) 

where 

f(k,x) = lk - 2c 2 h ^ X ~ X0) (3.12) 
and K c is the classical reflection factor of the model which is equal to 

R = ( {ik? + 2ifcyT + ^yT + ^ + 2(a + aQ \ / ik-2 \ 
c \{ik) 2 -2ifcvTT^v / TT^+2((To + a 1 ); \ik + 2j ' 
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The classical reflection factor (3.13) can be derived from the quantum reflection factor 
(2.8) by considering the classical limit i.e. when (3 — *■ 0. Because, in this limit || 
E = a + a\ and F Now following the idea introduced by Kim |15[] and 



developed by Corrigan \T%\, we may calculate the one loop quantum corrections to the 
classical reflection factor after perturbation calculation of the two-point function and 
then by finding the coefficient of e~' lk ( x + x ') j n the residue of the on-shell pole in the 
asymptotic region x,x' — > — oo. 

In order to calculate the one loop (0((3 2 )) quantum corrections to the classical 
reflection factor, we use the standard perturbation theory which is generalized 14-17 



to the affine Toda field theory on the half-line. In general, at (0(j3 2 )) there are three 
basic kinds of Feynman diagrams contribute to the two-point propagator of affine Toda 
field theory . These are shown in figure 1. Moreover, by inspection of the forms of the 
three point and four point couplings which we have found, it is clear that all types of 
these diagrams are involved in our problem. 






II 



III 



Figuer 1: Three basic Feynman diagrams in one loop order. 



In fact, when the boundary parameters are not equal then, the calculations corre- 
sponding to the one loop order in the sinh-Gordon model are lengthy and intricate. In 
the following sections we try to compute the contributions of types I and III diagrams 
to the reflection factor. The remaining diagrams will be treated elsewhere. Meanwhile, 
it is instructive to start with type III. 

4 Type III diagram (boundary-boundary) 

In this section we shall calculate the contribution of the type III diagram to the reflec- 
tion factor when both vertices are located at the boundary x — 0. We are led to the 
following integral 

— — (o"i cothx — cr tanhx ) 2 

x J J dtdt'G(x 1 ,t 1 ;0,t)G( y 0,t;0,t')G( y 0,t';0,t')G(0,t;x2,t 2 ) (4.1) 
in which the combinatorial factor has been taken into account. 
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Let us start by looking at the loop propagator G(0, t'; 0,f) which is equal to 
G(0,W) = if f%-T- n J , + • /(-^°) 

J J Z7T Z7T UJ ,Z — K ,z — 4 + ip 

x(j(V,Q) + K'(h>)f(-V,0j), (4.2) 

where 

IK + 2 

and K'(k') is the classical reflection factor ( |3.13| ). After some manipulation, we obtain 



G(0,t';0,t') 




duj' dk' 



2tt 2tt uj' 2 -k' 2 -A + ip 
2ik' (ik' -2 coth 2x ) 
X {ik') 2 - 2ik'y/TT^VTT^ + 2(a + aj ' ( ' ^ 



The above integral is clearly divergent however, the divergence can be removed by 
the infinite renormalization of the boundary term. In other words, considering the 
following relation 



2ik' (ik' - 2coth2x. 



oj 



{ik') 2 - 2ik'^TTa~ ^/TTa~ l + 2(a + ax) 

_ ^ (a/1 + <W1 + en - coth 2x ) - ((T + CTi) 

(ifc') 2 - 2ifc / vTT^vTT^ + 2(a + aj~' [ ' j 

it is seen that a minimal subtraction of the divergent part can be made by adding 
an appropriate counter term to the boundary, replace the logarithmically divergent 
integral by the finite part. Hence, 



G(0,t';0,lf) = 4i 




duj' dk' 1 



2tt 2tt uj' 2 -k' 2 -A + ip 
ik' (a/1 + cr \/l + <J\ - coth 2x ) - (a + a\) 



{ik') 2 - 2ik'^/TT^)V±T^~i + 2(a + a x ) 



(4.6) 



The integration over uj' may be performed by closing the contour into the upper half- 
plane and collecting a pole at uj' = \/k' 2 + 4 so that 

G(0 H- 0^=2/ dk ' 1 ^ (v /r +^v /I +^ - «)th2xo) - (gjg + gi) 
7 2tt a/F^+4 (2A; / ) 2 -2^ / v ^T^v / lT^+2(ao + a 1 ) 

In order to integrate over k', as before, one chooses a contour in the upper half-plane, 
however due to the branch cut the contour has to run around the cut line. Moreover 
we assume that the roots of the denominator of the integrand i.e. 2 cos ^ ao± ° jl ^ 7T are 
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positive, otherwise we may close the contour in the lower half-plane. Therefore, (|4.7|) 
converts to 

4 r !*y 1 V W 1 + <W1 + a i ~ coth 2a; o) + ( a o + a i) / 4 8 ^ 



and the above integral gives the following result 

G(0, £■ 0, t') = ~ cot a 7T - y cot aiTT. (4.9) 

Now it is convenient to calculate the time integral of the other middle propagator 
in ( |4.1|) which is equal to 



[ dt'G(0,t;0,t') = f f fdt' — -e-^ t -^ - 

J K ' ' ' ' J J J 2tt 2tt cu 2 - A; 2 - 4 + ip 

2ik (ik — 2 coth 2xn) 

x , , • (4.10) 

(z/e) 2 - 2iky/T+a^y/T+al + 2 ((To + o"i) 

Clearly, in the boundary-boundary contribution (|4. 1|) , it is seen that the if dependence 
is involved only in the propagator G(0,t;0,t') so, the integration over if gives us a 
Dirac delta function which means we must substitute zero for u and hence 

dk ( i \( 2ik (ik — 2 coth2xo) 



dt'G(0,t;0,t') 



2tt \-k 2 -AJ V. (ik) 2 - 2ik^l + a ^l + o\ + 2(a + o x ) 

(4.11) 

As we mentioned before, throughout this paper we assume that the roots of P(k) = 
(ik) 2 — 2iky/l + o" a/1 + o\ + 2(cr + cr\) which are equal to 2 cos ( ao± ^ n are positive, 
so the P(k) has no pole in the upper half-plane. Obviously, if the roots are negative 
then we can choose the contour in the lower half-plane in which no pole is inserted. 
Therefore, 

df G(0, f ; 0, t>) = , ^ 1 + mth2x °\ (4.12) 

and by substituting a = cosa vr and o\ = cosai7r, we obtain 

dt'G(0, t; 0, f) = -I — . (4.13) 

v ; 4cos^fcos^ v ; 

Up to now, the boundary-boundary contribution has the form 
i/3 2 (<Ti cothx — <Jq tanhx ) 2 (cto c °t a 7i + cti cot a^) 



32 cos *f- cos ^f- 



X 




dujt dk x f e -^i(*i-*) 
2tt 2tt u; 2 - fc 2 - 4 + z> 

+K 1 (A; 1 )/(-A; 1 ,x 1 )/(-A ;i ,0)e-^ 1 
, .du 2 dk 2 ie-^-^ f _ ik 

+K 2 (k 2 )f(-k 2 , x 2 )f(-k 2 , 0)e~^ ) . (4.14) 




First of all, it is necessary to perform the transformation k% — > —k\ in the first term 
of the first propagator. Secondly, integration over t ensures energy conservation at the 
interaction vertex and generates a Dirac delta function because of which we can set 
U\ = u>2- Moreover, it is better to define a new function as 

A(k, x) = f(-k, x)f(k, 0) + K(k)f(-k, x)f(-k, 0) (4.15) 

or, in an expanded form, 

A(k,x) 



ik + 2 coth 2xq ik + 2 coth 2(x — x ) 



+ 



ik + 2 ik - 2 

^)(zA; + 2co. , 



(ik + 2 cos (ao+ ° l)7r ) (ik + 2 cos (ao " ttl)7r s 



y ik - 2 cos (oo + ai)7r ) (ik - 2 cos 

ik — 2 coth 2xq ifc + 2 coth 2 (x — x ) 



zfc — 2 ik + 2 

then, the expression Q4.14| ) reduces to 



(4.16) 



i(3 (a i cothxo — <Jo tanhxo) (oo cot aovr + ai cot ai7r) 
32" cos 2f cos 




da^l dfei g-iwifo-ta) « * e -<fcia>i g-ifcai* 

27r 27r 27r - fef - 4 + ip u)\ - k\ - 4 + ip 



xA(ki,x 1 )A(k 2 ,x 2 ). (4.17) 

Obviously, what we need to do next is to integrate over the momenta k\ and k% and 
this may be achieved by closing the contours in the upper half-plane and considering 
the poles at k\ = k\ = k 2 = y/uf — 4. Note, the additional poles due to functions 
A(k\,xi) and A(k 2 ,x_) are not important because their contributions will be expo- 
nentially damped as x±,x 2 — > —00. Therefore, the boundary-boundary contribution 
is 

i(3 2 (a% coth x — <7o tanh x ) 2 ( a o cot aoii + a\ cot a\ii) 

cos ^f- cos ^f- 

g-i^Cti-ta) e -iki(x 1+ x 2 ) _____ A 0. u %i) A ^ x _y ( 41 g) 

Now recall the definition of the quantum reflection factor as the coefficient of 
e -ik(x+x') j n fog two-point Green function in the residue of the on-shell pole when 
x, x' — ► —00. Thus, the correction to the reflection factor from the type III (boundary- 
boundary) piece is 

iP 2 (<j_ coth^o — o"o tanhxo) 2 (a cot a vr + a\ cot aiiz) 
~1$2 cos *f- cos 

1 ( (i%_ + 2coth2x ) 2 nrjr/ * Aik 1 + 2co\h2x G )(ik l -2coth2x ) 
x—— — + 2K(ki)- 7—-^ 

2h \ (ikx + 2) 2 (ik t + 2) 2 

+ K 2 (k l ) ilh - 2 ^h2x ) 2 ' (4 ig) 

(ikt + 2) 2 



8 



5 Type III (boundary-bulk) 



This section deals with the determination of the contribution of the type III Feynman 
diagram to the classical reflection factor when one of the vertices corresponding to the 
loop is situated at the boundary and the other one is inside the bulk region. It is evident 
that in this case we have to take into account the bulk three point coupling C^ lk in 

(3) 

the corresponding vertex as well as the boundary three point coupling ClJ undary in the 
other vertex. Meanwhile, the combinatorial factor associated with the related Feynman 
diagram must be considered as a coefficient factor. Therefore, the contribution of the 
type III (boundary-bulk) to the reflection factor may be written as 

— 2/? 2 (o"i cothx — &o tanhx ) J J J dtdt'dxG(xi, t\] x, t)G(x, t; 0, t') 

xG(0, 0, t')G(x, t; x 2 , t 2 ) sinb^VWo) • (5.1) 

The propagator G(0, t'; 0, t') corresponding to the loop has been found in the previous 
section and is given by 

G(0, t'; 0, t') = -y cot a 7r - y cot ai7r. (5.2) 

The calculation of the other middle propagator i.e. G(x, t; 0, t') is the next step 
and clearly, the t' dependence in ( |5.1| ) is included only in this propagator. Hence, it is 
convenient to compute the following relation 

/* , G(»,f,o,o = /*'/ / g g e-M.-) ^ _ J_ 4 + , p ( m xW - k , „y» 

+A'(A)/(-<:,i)/(-A,0)e-' fa ).(5.3) 
Integrating over t' generates a Dirac delta function and so, 

fdt'G^W) = ij^ { _ k ^ A + ip) {mx)f(-k,0)e^ 

+K(k)f(-k,x)f(-k,0)e- ikx ) (5.4) 

and the residue theorem gives 

dt'G(x, t; 0, t') = — — (cq + ci coth 2(x — x )) , (5.5) 



where 



2 ( a o + ( 2 (a -ai)7r 



c = -ci = - ( 1 + tan^ 4 ) ( 1 + tan 2 ^— ; — — ) . (5.6) 



In order to check the above result, if we set x = in ( |5.5|) then it will be equal to 
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Up to now the type III (boundary-bulk) contribution take the following form, of 
course, after integrating over t: 

(3 2 co((Ji coth xq — do tanh x ) (oo cot a ir + a\ cot aiir) 

* L * / / 1 £ ^""^ ^W^TTP (/(*.• ,)e*<"' 

+K 1 (k 1 )f(-k 1 ,x 1 )f(-k 1 , x)e~^ x+ ^ 

(1 - coth 2(x - x )) sinh(v / 2/?0 o 



+K 2 (k 2 )f(-k 2 ,x)f(-k 2 ,x 2 )e~ <b ^). (5.7) 



By multiplying the two propagators in ( |5.7|) by each other, it is clear that one obtains 
four pole pieces and, as far as the integration over x is concerned, if we do the integration 
over x on one of them then, obviously the other three pole pieces could be done in the 
same manner. Hence, in what follows it is sufficient to treat only one of them and, 
meanwhile, keeping those terms which are functions of x, we are led to the following 
complicated integral 



/dx exp {2 + i(k 2 — ki)x}f(—ki, x)f(k 2 , x) sinh(v^2/30o) ( 1 — coth 2( 
-co ^ 



X - X 0/ 



(5.* 

After some substitutions and collecting together powers of coth2(x — xq) we obtain 



/dx exp {2 + i(k 2 - ki)x} sinh(V2/30 o ) ( ~hh 



(iki-2)(ik 2 + 2) 

+ (2ik 2 - 2iki + k x k 2 ) coth 2(x - x ) + (2ik x - 2ik 2 - 4) coth 2 2(x - x ) 

+4coth 3 2(a;-xo) N ). (5.9) 



It is clear that in order to solve the above integral, it is necessary to manipulate the 
following integrals 

/ dx exp {2 + i(k 2 -k 1 )x} sinh(V2/50 o ) coth" 2(x-x Q ), (5.10) 

J — oo 

where, n — 0, 1, 2, 3. 

In fact in Appendix A, we have found the integrals ( |5.10| ) and the solutions of them 
are expressed in terms of hypergeometric functions. So, using the formulae in Appendix 
A and simplifying, we find that J in (5.9) can be rewritten 



F{ki, k 2 ) 



3hk 2 + Aik 2 - Aikx + 8 1 



3k\k 2 + Qik 2 — Qiki + 13 cosh2xo 



sinh 2xq 
10 
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sinh 2x, 



o 



ik 2 — ik\ + 2 cosh 2 2x + 1 1 cosh 3 2x + 5 cosh 2x 



o 



3 sinh 3 2x 6 sinh 2x 

\2ik x k 2 - 16k 2 + I6h + 40i - (fc 2 - fci) (9fcifc2 + 16ife 2 - 16ifci + 34) 

3(fc 2 - fci - if) 

x e~ 2x °F (l, % -{k 2 - h) + 1, % -{k 2 - h) + 2, e" 4 ^ 

\2ik x k 2 - 48fc 2 + 48fc x + 136z - (k 2 - fci)(6fcifc 2 + 28ik 2 - 28ih + 84) 

3(fc 2 - fci - 8i) 

x e- fo0 F ^2, ^(fc 2 - fcx) + 2, i(fc 2 - fc x ) + 3, e" 4 ^ 

32fc 2 - 32*4 - 192i + (fc 2 - fci)(16ifc 2 - 16zfci + 104) 
+ 3(fc 2 - fci - 12i) 



- 10x °F ( 3, i(fc 2 -fc 1 )+ 3, 1(^-^)4- 4, 



^16fc 2 - 16fci - 32i 
(fc 2 — fci — 16i) 

x e~ 14:Eo F (4, ^(fc 2 - fci) + 4, ^(fc 2 - fci) + 5, e" 4 ^ . (5.11) 

Now regarding (fT7), after doing the transformation fci — > —fci in the first term of the 
first propagator, all that remains is to integrate over the momenta fci and fc 2 and this 
can be achieved by closing the contours in the upper half-plane and considering poles 
at fci = fci = fc 2 = \f oj\ — 4. The extra poles in the four functions .F(±fci, ±fc 2 ) are not 
important because their contributions will be discounted when x% and x 2 go to — oo. 
Let us write down the type III (boundary-bulk) contribution to the reflection factor 



2 ton Oo+ai)-"- fQT1 (ao-ai)ir 

(do cot aoTT + cli cot a\7r) 



(3 2 tan tan j 

„„„ Q17T 



COS ^f- COS ■ 2 

g-iwi^-fa) e -tfc!(xi+x a ) 1 + 2coth2(xi -x ) ik 1 + 2coth2(x 2 -x Q ) 

2vr ' (2fci) 2 ifci - 2 ffci - 2 

x <j k) ~ —~ —« ^i(fci)^(fci, fci) 

\zfci + 2) 2 (zfci + 2)(zfci-2) 

-71 ^71 -KiChm-h, -fci) + ri 1 — Kfikmk -fci) 1 (5.12) 

[iki + 2){ik\ — 2) [iki—2y \ 

Now looking at the function !F(ki, fc 2 ) given by ( |5.11| ), let us show the detailed forms 
of F{— k±, fci), ^"(fci, fci),J r (— fci, —fci) and .F(fci, —fci). In fact, 

r , t f N 3fc? + 8 1 3fc? + 13 cosh2x 

J 7 fci, fci = 5 ■ r + 7 . u2o 

3 smh2x 6 smh 2x 

2 cosh 2 2x + 1 1 cosh 3 2x + 5 cosh 2x 

3 sinh 3 2xo 6 sinh 4 2x 
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+ 3fe i + 10 e -2so F (i >;L) 2, e- 4 *°) 

+ 3fc l+ 34 e ^o F(2)2i3)e -4xo ) 

6 

16 

+ — e~ l0xo F(3,3A, e~ Axo ) 
3 

+2e~ Ux ° F(4,4,5,e" 4a;o ). (5.13) 

It can be easily verified that 

F(-k u -k)=r{k,ki), (5-14) 

and 

w f f , 3fc? — 8iA;i — 8 1 3A; 2 - 12ijfei - 13 cosh 2x 

Si—kijKi) = : 1 5 

3 sinh 2x 6 sinh 2x 

2ik\ + 2 cosh 2 2x 1 cosh 3 2xq + 5 cosh 2xq 



3 sinh 2x 6 sinh 2x 

9kf - 38ikj - 50h + 20i 



e-^F(l,-k 1 + l,-fc 1 + 2,e- te «») 
3(fci — 2z) ^ ^ 

6fc 3 - 62^ 2 - 132k + 68. ^ i , i £ _ 



e- ta0 F(2,-fci + 2,-^ + 3, e" 



3(/fc! - 4z) 2 2 

32ifc 2 + 136*1 - 96i _ Wx i f i 



3(fci - 6i) 2 2 ; 

+ 16^-16i e _ 14 , o i ^ + ^ t_ h + ^ e _ 4xoy (5 15) 

(hi — 8i) 2 2 

Finally T{ki, —ki) can be obtained from T{—ki, k\) after setting fcj — > —k±. 

6 Type III (bulk-boundary) 

In this section we study the quantum correction to the classical reflection factor due 
to the contribution of the type III Feynman diagram, when the vertex associated with 
the loop is located at the bulk region and the other vertex coincides with the boundary. 
The associated contribution is given by 

C = — 20 1 {o\ cothx — cr tanhx ) J J J dtdt'dx'G(xi, ti, 0, t)G(0, t; x', t') 

xG(x', t'; x', t')G{0, t; x 2 , t 2 ) sinh(V2/30 o )- (6.1) 

The following relation which is some part of the contribution ( |6.1|) , can be derived 
independently from the remaining part 

Ci = J dtG{x 1 ,t 1 ;0,t)G(0,t;x 2 ,t 2 ) (6.2) 
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or 



Cl = /dt //| 1 ^ e -^,-0_L_( /( , 1 , Il)/( _, 1 ,o ) e-. 




+K 1 (k 1 )f(-k 1 ,x 1 )f(-k 1 ,0)e- ik ^ 

x / Ippe-**™ 2 p^— - (f(k 2 ,0)f(-k 2 ,x 2 )e^ 
2tt 2tt — K2 — 4 + ip \ 

-K 2 (k 2 )f(-k 2 , 0)f(-k 2 , x 2 )e-' lk2X ^ . (6.3) 




First of all, it is necessary to set ki — > —k\ in the first term of the first propagator. 
Secondly, integration over t leads to the substitution of u 2 = Finally integration 
over the momenta k\ and k 2 , as before, may be done immediately by closing the 
contour in the upper half-plane and looking at the poles at k\ = k\ = k 2 = \J u 2 — A 
and ignoring all the other poles as their contributions vanish rapidly as xi,X2 — > —00. 
Therefore, 

_± g-iwiCti-ta) g-tfci^+M) __ A ^ h Xl )A(k u x 2 ), (6.4) 

where 

A(h, Xl ) = fi-k^x^fCh, 0) + K(k)f(-h, xi)f(—ki, 0). (6.5) 

So, our next problem is to calculate the following integral which is the remaining part 
of the contribution. 

dt'dx'G(0, t- x', t')G(x\ t'; x', f) sinh(V2/30 o ). (6.6) 




Obviously, this part will be appeared as a constant and it must be multiplied by 
Clearly, the time variable t' appears only in one of the propagator i.e. in G(0,t\ x',t'). 
On the other hand, this propagator along with integration over if has been obtained in 
the previous section. Hence, 

J dt'G(0, t; x', £') = ^e 2x 'c (1 - coth 2{x' - x Q )) , (6.7) 



where 

c = - ( 1 + tan 2 v ^" u ' ~ u " ) ( 1 + tan 2 / 17 " ) . (6.{ 



2 (a + ai)7i\ ( 2 (a — ai)ir 



Therefore, (|6.6|) reduces to 

Y J dx'e 2x ' smh(V2(3(f) )(l-coth2(x' - x ))G(x',t';x',t'), (6.9) 



where 



GW.t':.r'.t>) = ! I J [fik'.s')fi-k'.s' 

2it 2it uj z — k — 4 + ip 




+K(k')(f(-k',x')(f(-k',x>)e- 2lk ' x ') (6.10) 
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or after integration over u' 



1 f dV 1 

G(x', t'; x', t') = - — - ( f(k', x')f(-k', x') 

2 J 2n ^k' 2 + 4 1 



+iT(A;')(/(-A; / ,x / )(/(-fc , ,x')e- MV ). (6.11) 



In fact, the above integrand has two parts, the first part can be easily manipulated but 
the other part which includes the exponential term is hard to calculate and we prefer 
to leave the computation of that part for later. Let us look at the first part of the loop 
propagator. The integral of this part is logarithmically divergent. Nevertheless, this 
divergence can be removed by an infinite renormalization of the mass parameter in the 
bulk potential. Then, doing the integration over k', we obtain 

i f^m-^^- ^-^y-^ . (6,2, 



2 J 2vr ^/WTl ' ' 2tt 

To sum up, the integral ( p.9| ) reduces to 

r-0 



16n 



/ dx' e 2x ' ( 1 - coth 2 2(x' - x Q ) ) ( 1 - coth 2( x' — Xo )) sinh(\/2^ ) 

/0 p i j f -i 
dx' / ( 1 - coth2(x' - xo)) sinh(VWo) e 2x ' 

.oo y 27r VA;' 2 + 4 v / 



16 

{ik'f + 2ifcVl + <W1 + ^1 + 2(ct + o-i) (i# + 2coth2(x' - x )) _ 2;7 , v 



(zA;') 2 - 2zA; , v ^Ta " v / lT^ + 2(a + <n) (i# + 2)(ik' - 2) 



e- 2lk x K6.13) 



The above relation has two parts and the first part which is a single integral can be 
performed by means of the formulae in Appendix A and we write down only the solution 
of this part which is expressed in terms of hypergeometric functions, that is, 

C 2 = / dx'e 2x ' (1 -coth 2 2{x - x )) (1 -coth 2{x - x )) sinh(V2/30 o ) 

167T J_ 00 

1 ^l + tan 2 ^±^Vl+tan 2(a °- ai)7r 



16tt V 4 y V 4 

1 cosh 2xq cosh 2 2xq + 1 cosh 3 2xo + 5 cosh 2xc 



x 



3sinh2x 24 sinh 2x 6sinh J 2x 24 sinh 2x 

--e~ 2xo F {1,1,2, e'^ ) 

+ H e -6-o F ( 2) 2,3, e - 4:Co ) 

+^-e- 10xo F (3,3,4, e~ 4xo ) 

+^e- Uxo F '(4,4,5, e" 4a:o )| . (6.14) 
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So, in connection with the type III (bulk-boundary) contribution, the remaining inte- 
gral is 

Cs = S/° dx ' I % 7W+t ( x ~ coth2(x ' ~ Xo) ) sinh ^^ e2x ' 

(ik' + 2 cos (oo+ 9 ai)7r > ) ( ijf + 2 cos (a °~ ai)7r N ) f ifc' + 2 coth 2(x' - x )) 
X V 2 / V 2 / V V c -2tfcv ( 615 ) 

^'-2cos^^ViA;'-2cos^i)£^ (i# + 2)(i# - 2) 

As we mentioned before, it is more convenient to integrate over x' then afterwards over 
k! . Since to integrate over k' first is a difficult problem. Let us do partial fraction de- 
composition for the rational function in ( |6 . 1 5[ ) . Obviously we will have four elementary 
partial fraction including 

1 111 



(oo+oi)7r ^ ' ( -y _ 9 (ao-oi>7iA ' ik' + 2 ' ik' — 2 



ik> - 2 cos w+£i£L j ^jfe/ _ 2 cos ■ , 

Now, in what follows we perform the calculations in detail for one of them, for 
example, j— / 1 (ao+ai)7r x due to the fact that for all the others the computations are 

I zfc' — 2 cos 2 ) 

similar except that cos ( - ao+ ^ 1T i s replaced by one of cos ("o"" 1 ) 7 ^ _i ; i ; respectively. 
What we need to do is to calculate the following : 

j- aT1 2 (ao+ai)7r _ ,2 (ao+Q]>r 

ld,u 4 tOL 4 / a vr ai7r 

cot cot 



cos 2 (ao+ 4 ai) " cos 2 (a °~ ai) " '22 



if f dk' 1 / \ 
x — / dx' / ( 1 — coth 2(x' — xp) ) sinh(V2d(b ) e 



(2-2iJfe')a:' 



(a + ai)7r\ 2 / 1 



x coth 2(x'-x ) + cos ^ ^- , - , . (6.16) 

The integration over x' may be done by using the formulae in Appendix A and gives 
I dx' e {2 - 2ik ' )x ' (l - coth 2(x' - x )) sinh(V2/30 o ) 

J — oo 

x (coth 2(x'-x )+ cos ( a ° + fll)7r ) 2 

= L(a ,Oi) 

+ E { ^^r e " (2+4(n " 1))ao F ^ - ^ + n » -- 2 k ' +n+l > e_ "°) } ^ 

where 



2 (a + cli)ii 4 (ao + ai)7r 2 



L(a , ai) = cos cos 



2 3 2 3 / sinh 2x 
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2 (ao + a,i)ir (ao + aiV 13\ cosh2£ 

+ l 2 C0S ^ C ° S 2 + 24J iliTi?^ 

1 (ao + ai)n 1 \ cosh 2 2xq + 1 

— I - COS 5 

3 2 QJ sinh 3 2x 

1 cosh 3 2xq + 5 cosh 2x , n , oN 

+ 1 (6.18) 

24 sinh 4 2a;o 

and the coefficients A n , B n , n = 1,2,3,4 are constants which in fact only depend on 
cog (aoj^iK _ n ow the final calculation is to integrate over k' and it is evident that in 
order to do this, we have to convert the hypergeometric functions to infinite series. 
Considering the equation ( |A.9| ) in Appendix A, we conclude that 



F% -\k' + 1, -\k' + 2, e- 4 *°) _ ~ e- 4 ^° 

k' + 2i ^ k' + i(2 + 2n)' ^ ' ' 

If we differentiate both sides of the above relation with respect to xq, then the following 
identity may be derived 



F(2,-\k' + 2, -{k' + 3,e" 



ix \ JO. np- 4 ( n - 1 ) a: 



fc' + 4z ^k' + i(2 + 2n 

n=l 

In the same way, one obtains 

F(3, -ifc' + 3,-ifc / + 4, e- 4 * ) _ 1 n(n - i) e - 4 ("" 2 )^ 

/c' + 6i ~2!^ fc' + z(2 + 2n) ^' ^ 

n=2 v y 

and 

F(4, -\k' + 4, -ifc 7 + 5, e~ 4x °) _ 1 A n(n- l)(n- 2 )e- 4 ("- 3 ) a:o 

A;' + 8* ~3!^ fc' + i(2 + 2n) ' ^' ^ 

n=3 ' 

Now if we substitute (6.19), (6.20), (6.21) and (6.22) into (6.17), all that remains 
in connection with the contribution (6.16) is the integration over k' . Obviously we 
encounter integrals of the form 



dk' 1 \ f Ak' + B 



-ooVWTl \ik> -2cos^il?E/ \k' + i(2 + 2n) 



(6.23) 



and the k' integration may be performed by closing the contour in the upper half-plane 
and onto the branch cut which stretches from k' = 2i to infinity along the imaginary 
axis. In fact, leaving the integrals along the branch cut to be evaluated later, we obtain 
the required formula 

dk' ( I \ / Ak' + B 



VW+4 \ik'- 2 cos (a °+ ai) " / \k' + i(2 + 2n) 
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2 cos (ao + ai) M + iB 



(oo+oi)7r 
2 



2n + 2 - 2 cos 



(ao+ai)7r 



sin 



((2n + 2)A + iB) 



(ao+ai)ir 
2 

1 



In 



n 



+ 1 - Vn 2 + 2n 



^2n + 2 - 2 cos (a °+ ai)7r j 2\Jn 2 + 2n | n + 1 + Vn 2 + 2n 



. (6.24) 



Note (|6.24j) is valid when n ^ 0, on the other hand if n = then one may find 



dk' 



\ik'- 2 cos (Q °+ ai) " 
2cos (ao+ 2 ai) M + iB| i^,^... 



fc' + 2i 

(2A + z\B) 



2 - 2 cos 



(ao+ai)7r 



(ao+ai)7r 



sm 



2 - 2 cos 



(ao+ai)7r 



(6.25) 



Now we are in a position to write down (|6.16|) or, in fact, the integral (6.15) 



t ari 2 ^llE _ C ot 2 



32vT cos 2 (°0+°iK CQS 2 (°0-ai)ir 



cot — — cot — — 



(ao+ai)7r 
2 



sm 
1 

+ I2 e 



(oo+ai) 



— L(a , ai, 



2 

-2xo 



2Ax + zSi 2 cos I^ilEAi + iBj ^±^L 



2 - 2 cos (a °+ ai)7r 



2 - 2 cos (a °+ ai)7r 



(ao+ai)7r 



sm 



3+ai)7r 



-(2+4n)z 



..jsm^ ( 2 n + 2-2cos^±^ 



a + aivr 

2 COS A\ + 2±5l 



n 



+ Yj I 2 cos 



(a + ai)vr 



A 2 + iS 2 + 



n(n — 1) 
2! 

n(n — l)(n — 2) 



3! 



(an + aiW , 
2cos l 2 1J A 3 + iB 3 

(an + GhW . „ 
2 cos v " ; A 4 + iB 4 



e - (2+4n)rr 



t (2n + 2-2cos (a ° + 2 ai)7r ) 2 v^ 2 + 2n 
(2n + 2)Ai + iSi 



In 

n 
1! 



n 



+ 1 - V^ 2 + 2n 



n 



+ 1 + Vn 2 + 2n 



{2n + 2)A 2 + i£ 2 



n(n — 1) 
2! 



(2n + 2)A 3 + iBz 



+ 



n(n — l)(n — 2) 
3! 



(2n + 2)A A + zfi 4 



+ other pole pieces. (6.26) 

Note, in the above expression all the series are convergent. As we mentioned before, 
( |6.26| ) must be considered (after adding to (6.14)) as a coefficient factor of ( |6.4| ) in 
order to constitute the type III (bulk-boundary) contribution i.e.: 

C = d (C 2 + C 3 ) . (6.27) 
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7 Type I diagram 



In this section we calculate the contribution of the type I Feynman diagram to the 
classical reflection factor when the vertex is placed inside the bulk region. In fact, 
when the vertex is located at the boundary then, the corresponding contribution has 



been found |24| and is given by 

if3 2 

— — (<7i coth xo + o"o tanhxo)(ao cot aovr + a% cot axiz) 
8 

x / — e -*"(*i-ta) e -&(*i+*2) ik + 2coth2 ( x i -xp) ifc + 2coth2(x 2 -xq) 
J 2tt ' p{k) p(k) 

Clearly, in our case the bulk four point coupling should be considered in the in- 
teraction vertex. Moreover, as before, the combinatorial factor associated with this 
diagram will appear as a coefficient. Hence, the contribution has the form 



/oo fO 

dt / dxG(xi,ti, x,t)G(x,t; x,t)G(x,t; X2,t2) cosh(V2f3(f) ), (7.2) 
■oo J — oo 

where 

cosh(v / 2/30 o ) = (2 coth 2 2(x - x ) - l) . (7.3) 
In the previous section, we simplified the loop propagator to 

(1 - coth 2 2(x - x )) 
G(x,t;x,t) = -A ^ °JL 

1 / ok_ K(k")f(-k", x)f(-k", x) e- 2lk " x . (7.4) 

2 J 2tt y/k" 2 + 4 

Also the integral part of the loop Green function is hard enough to evaluate and we 
found out that it is better to do this integration during the final stage. Now, let us 
rewrite the contribution (|7.2|) in the expanded form 



o 




duj dk 



-Ai(3 2 dt dx e -M*i-*) (f(k,x 1 )f(-k,x)e ik ^ 

27T 27T uj — k — 4 + ?p V 



+K(k)f(-k, x x )f(-k, x)e- ik{xi+x) ) cosh(VW, 



(1 - coth 2 2(x - x )) , 1 fdk" 1 ,„ w 



x < -- — 2V — - + 2 y ar^Ti^(-^^-^^ 

>< / / g g ^ - * - 4 + *p ( «>^' 

+K'(k , )f(-k', x)f(-k', x 2 )e~ M ^) . (7.5) 

Looking at ( [7. 5[ ) , one can predict that the calculations will be lengthy and intri- 
cate. The starting point is to do the t integration which allows the substitution uj = u'. 
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Secondly, it is necessary to perform a transformation k — > —k in the first term of the 
first propagator. Moreover, if we multiply the first and the third propagator with each 
other, then obviously we will have four pole pieces and fortunately if we do the calcu- 
lation for one of them (for example the first one), then the calculations corresponding 
to the other three pole pieces may be treated similarly with k + k' replaced by one of 
k — k', —k + k' and —k — k! . Because of this in what follows we follow the problem 
only for one pole piece. Hence our problem is, in fact, the following integral 



V = -^ f° dxf t d ^ d ± e -^) 
J-oo J J 2tt 2tt 



to 2 — k 2 — 4 + ip 
x cosh(v / 2 / 90 o )/(-A;, xi)f(k, x) e - ik(xi - x) 

(l-coth 2 2(a;-a;o)) , 1 f dk" 1 WLffw , u , ^^_ 2ik „ x 



x ^ —+~ / -7. K(k")f(-k",x)f(-k",x)e 

1 2tt 2 J 2tt ^k" 2 + 4 



x 



In fact, the above contribution has two parts. The first part, in which the integral of 
the middle momentum (k") is not involved, can be calculated by means of the formulae 
in Appendix B and we call this part T>\. Let us write down the solution of this part. 
This contribution is expressed in terms of the hypergeometric function as : 



Vl = !^y*^ e -Mti-te) e -^i+- 2 ) /( _^ Xi)/( _^ X2) j_ 

F^2,^ + l,^ + 2,e- 4 ^ 



n 2 

Ai-Ak- ik 2 
k-2i 



x < ; e " 4xo 



A8i - 40k - 8ik 2 8r ! ? n i t 

+ s e~ 8x ° F 3, -k + 2, -k + 3, e~ 4x ° 

k - Ai V 2 2 

176z — 96A; — 8ik 2 19 „ „ / i f i f . 

+ s e" 12xo F 4, -jfc + 3, -k + 4, e" 4:ro 

k — 6i V 2 2 

+ 25y-64fc e _ lfaoF / j fe + 4 j fe + 5>e _ teo 



fc - 8i V 2 2 



Now it is better for the second part, which we call T> 2 , to integrate first over x 
then over k" . Meanwhile, before starting the integration, it is useful to note that if we 
do the partial fraction decomposition for K"(k")f(—k",x)f(—k",x), then we will have 
four elementary partial fractions as 



(ijfc//_2cos*22±p*L)' (ik" -2 cos lk " + r lk "~ 2 ' 
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As before, in the remaining section we continue the computations in detail for one of 
them (for example, 1 (a[)+ai) ^x ) because the calculations corresponding to the 

I zfc" — 2 cos 2 I 

other three elementary partial fractions can be done in the same manner just by the 
substitution of cos ( ao + ai ^ >7r by one Q f cos l - a °~^' n ^ _i ) i ; respectively. So, our problem 
reduces to this integral 

2*/? 2 cot ^ cot ^ f tan 2 ^±^1 _ co t 2 ^ ± ^ 



2 2 



o 




duj dk „•, ,,, , s i 



x dx e -M*i-fa) 



2tt 2tt lu 2 - k 2 - 4 + ip 

x cosh(VWo)/(-fc, xx)f(k, x) e ik{xi ~ x) 

dk" 1 e~ 2ik "* ( (ao + a 1 )7^^ 2 

x / ~oZ~ /ttto : r ~T7 ^ coth2(x-x ) + cos 



X 



2tt VA;" 2 + 4 ijfe" - 2 cos (a °+ ai)7r V " ' ' ' 2 

^- 2 , /2 , , ■ * 2 ) e^"- ). (7. 

27T CtT — « w — 4 + ?p 



The integration over x can be evaluated by means of the formulae given in Appendix 
B as 



f° dx el (k +k >-2k"y, cosh (72/30 o )/(A;, x)f(k', x) Loth. 2{x - x ) + cos (a ° + Ql 

J-oc V 2 

1 (A' Q kk' + B' {k + fc') + C' ) 



2 

)7T ' 



(ik + 2){ik' + 2) (k + k' — 2k") 



1 / (A' n kk' + -B^(& + k') + C^j) e _4( n _i) x . 



(ik + 2)(ik' + 2) ^ \ (k + k' - 2k" - 4(rc - 1) 

xF^(Hi'- 2Jfe") + n — l,^(fe + A:' - 2fc") + n, e" 4 ^ J , (7.9) 

where the coefficients A' n , B' n , C' n ; n = 0,1,. .,6 are constants and depend only on 
cos ( a °+ ai ) 7r , i n fact A' 5 , A' 6 and B' 6 are zero. Now the subsequent calculation is to 
integrate over k" and it is clear that to do this, it is necessary to convert the hyperge- 
ometric function to an infinite series. Looking at (|B.9|), we may write down 



F(l, \(k + k'- 2k"), \{k + k'- 2k") + 1, e- 4 *°) " e~ 4 -° 

k + k 1 - 2k" ^ k + k 1 - 2k" - Ani ^ ' ' 

n=0 

and by differentiating both sides of the above relation with respect to Xq, then we 
obtain 

F(2, {{k + k' - 2k") + 1, \{k + k' - 2k") + 2, e" 4 * ) g ne -±(n-i) xo 



k + k' - 2k" -U ^ k + k' - 2k" - Ani 

n=l 



(7.11) 
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Similarly one can derive the infinite series forms of the other hypergeometric functions 

as 

F(3, {(k + k'- 2k") + 2, {(k + k'- 2k") + 3, e~ 4x ») _ 1 ~ njn - i) e -^-2)x 
k + k' - 2k" -8i ~ 2! ^ k + k' - 2k" - 4ni ' 

n=2 

(7.12) 

F(4, £(*; + #- 2ife") + 3, ±(Jfe + k' - 2k") + 4, e - 4:C0 ) 



k + k 1 - 2k" - 12i 
1 ^ n(n - 1) (n - 2)e- 4 ( n " 3 ) :Eo 
~3!^ k + k' - 2k" - Ani ' ( ' ^ 

71=3 



and 



F(5, {{k + k' - 2k") + 4, {(k + k' - 2k") + 5, e~ 4x °) 
k + k' - 2k" - m 

1 ^ n(n - l)(ra - 2)(ra - 3) e - 4 ( n - 4 ^ 
" 4! ^ k + k' - 2k" - 4m 

n=4 



F(6, + fc' - 2k") + 5, i(ib + k' - 2k") + 6, e~ 4x °) 



(7.14) 



k + k' - 2k" - 20i 

1 y> ti(ti - l)(w - 2)(w - 3)(n - 4)e- 4(ra ~ 5):ro 
~5!^ k + k! - 2k" - 4m' ' ^' ^ 

n=5 

Let us substitute (7.10), (7.11), (7.12), (7.13), (7.14) and (7.15) in (7.9) and obvi- 
ously what remains in connection with the contribution (7.8), are the integrations over 
k", k' and k . As before in previous sections, in order to integrate over the momenta k 
and k', it is sufficient to close the contours in the upper half-plane and pick up poles 
at k — k — k' — \fuj 2 — 4 as all the other poles' contributions will be exponentially 
damped when x, x' — > — oo . Meanwhile, the integration over k" is of the form 



/ 



dk " 1 1 (7.16) 



V^^Tl (i k » _ 2 cos (k + k' — 2k" - 4m) 



Now to manipulate the integral, let us choose the contour in the upper half-plane, 
taking care of the branch cut which runs from k" = 2i to infinity along the imaginary 
axis. Clearly this integral reduces to the integrals along the branch cut and we obtain 

dk" 1 1 



Vl^TA - 2 cos (ao+ 2 ai)7r ) (k + k'- 2k" - 4m) 



(oo+oi)7r 
2 



(k + k' + U cos (a °+ ai)7r - 4m) V sin 



(ao+ai)n 
2 
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2% 

i fc ±^ + 4-4n 2 -2m(fc + A; / ) 



i -^ ±1 + n + fV^T^ + 4 - 4n 2 - 2ni(k + k>) 



1 + + n - iA/^r^ + 4 - 4n 2 - 2m'(A; + fc') 



When n = 0, the above formula is simplified much more, especially after doing the 
integration over /c and fc' and using the fact that k = k = k' = 2 sinh 9 . So the 
following formula can be obtained 

dk" 1 1 



Vk" 2 + 4 ^// _ 2 cos is^^j ^ 2 jfc - 2A;" 



(ao+ai)7r „ 

+ -«— (7.i8) 



Now, (7.8) or more generally the contribution T> 2 can be written as: 

i(3 2 a 7r a x i\ ( 2 (a + ai)vr 2 (a + a x )i\ 

V 2 = - — cot cot tan cot 

4tt 2 2 V 4 4 




E 



e 



-4nx / ( a o+ a i ) 7r 



2 



+ 2) 2 ^ ( 2 * + At cos - 4m) \ sin 



2i , f 1 + | fc + n + | a/A; 2 + 4 - 4n 2 - 4m fc 



a/A; 2 + 4 - 4n 2 - 4m ll + |A; + n- |\/A; 2 + 4- 4n 2 - 4mfc 



2 1 2 

(^A; 2 + 2B[k + CO + n(A' 2 P + 2B' 2 k + C£) + ~ ^ (A' 3 P + 2fl£jfe + C£) 
+ n(n-l)(n-2) ( ^ + ^ + ^ + n(n - l)(n- 2)(n - 3) ^ + ^ 

n(n - l)(n - 2)(n - 3)(ra - 4) 
+ 5! ° 6 

+ other pole pieces. (7-19) 

Firstly, in order to check the above solution, if we set and consider the other pole 

pieces then, we can derive the formula (3.10) in reference |p^| . As we mentioned before, 
the calculation of this reference is based on the case when the boundary parameters are 
equal. Secondly, in this solution, we verified that the term which depends explicitly on 
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the rapidity of the particle (9) is cancelled by counterpart terms in the other pole pieces. 
It is evident that if we add the expressions (7.7) and (7.19) then, the contribution (7.6) 
will be obtained i.e. V — V x + T> 2 . 



8 Discussion 

Affine Toda field theory on the whole line is an exactly solvable theory for which 
the S-matrices have been formulated. However, when a boundary is present then the 
boundary S-matrices of the theory i.e. the reflection factors, have not been completely 
found. The bootstrap technique does not uniquely determine the reflection factors. 
Fortunately perturbation theory provides the link between the expressions for the re- 
flection factors which come from the bootstrap equations and the boundary parameters. 
Nevertheless, this method normally involves complicated calculations. 

In this paper the quantum reflection factor for the affine Toda field theory or 
sinh-Gordon model with integrable boundary conditions has been studied in low order 
perturbation theory when o"o 7^ 0i- It is found that at one loop order the quantum 
corrections to the classical reflection factor of the model can be expressed in terms of 
hypergeometric functions for most of the related Feynman diagrams. Although there 
is still some work to do to calculate the contributions of the remaining diagrams, it 
is understood that the provided procedure and some formalisms may be followed for 
them. 

The calculations corresponding to the type II Feynman diagram which are not 
carried out in this paper, are more difficult than the others. In this case the two 
middle propagators are exactly the same and this fact influences the difficulty of the 
computations. However some formulae that have been presented here, could be helpful 
for the remaining diagram. For example, consider the contribution of the type II 
(boundary-bulk) diagram: 

— 2f3 2 {a\ coth^o — (To tanh^o) J J J dtdt'dxG(xi,ti;x,t)G(x,t;0,t') 

xG(x, t; 0, t')G(0, t'; x 2 , t 2 ) smh(V2~P<j> ). (8.1) 
Now as far as the integration over x is concerned we should obtain the following integrals 

/ dxe i{k+k '' kl)x smh(V2(3(j)o)coth n 2(x -x ), (8.2) 

J — oo 

where n=0,l,2,3. It is better to solve: 

r° 

/ dx exp {r + i(k + k' — h)x} sinh("\/Wo) coth n 2(x - x ) (8.3) 

J — oo 

in which r is a small positive quantity and will be taken to zero later. In fact, the 
relation (8.3) is very similar to the formula (A.l) in Appendix A. So, following the 
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same procedure that have been followed in Appendix A, one can find the solution of 
(8.2) when n = as: 



-oo 



1 



sinh 2xq 

2(k + k'-k 1 ) 
k + k' - h - 2% 

x e~ 2 *°F(l, l(k + k' — h) + l, l -(k + k'- h) + |, e" 4 ^). (8.4) 

Then, the solutions of (8.2) for n = 1, 2, 3 can be derived exactly in according to the 
Appendix A terms. But, this is not all of the problem. As we mentioned before, in 
type II diagram double Green functions cause the middle momenta to be linked to each 
other in a complicated way and the calculations become more intricate. Actually this 
diagram must be studied in three cases depending on the interaction vertices being 
located in the bulk region or at the boundary. Moreover because of the symmetry, 
the contribution of the type II (boundary-bulk) diagram is the same as the type II 
(bulk-boundary) one. 

When the boundary parameters are equal then, only the type I diagram is involved 



in the theory. As we mentioned before, in this special case [|T4]] the quantum corrections 
to the classical reflection of the model have been found and Ghoshal's formula for the 
lightest breather is checked perturbatively to 0{(3 2 ). In our case, we realised that the 
contribution of the type I (bulk) reduces to the special case. Taking (7.7) and (7.19) 
expressions, if we put er = o\ then, we obtain the same result as reference [14]] and 
this is a check on our calculations. Moreover, when cr ^ <J\ the following expressions 
for E and F in Ghoshal's formula (2.8) have been conjectured [21J to be: 

E = (a + a 1 )(l- B/2) F = (a — ai)(l — B/2). (8.5) 

So, it will be interesting to check the above conjecture after finding the contributions of 
the remaining diagrams and adding the results all together. This will lead to a deeper 
understanding of the quantum integrability of the theory. However, it is necessary to 
find simplifications of the contributions when they add among themselves in order to 
get Ghoshal's formula. 
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Appendix A 

In this Appendix we obtain such integrals 

f° 

S n = dxe {2+lk)x smh{V2(3(i)v)cot\i n 2{x-xv) (A.l) 

J — oo 

in which n — 0, 1, 2, 3, O is the background solution to the equation of field so that 
sinh(v / 2/50o) is proportional to the bulk three point coupling which is given by 

sinh( VWo) = 2 cosh 2(x - x Q ) (coth 2 2 (x - x ) - l) . (A.2) 

Let us start with the simplest case when n = 0. Using ( |A.2| ), we have 



o / 1 



So = ~ — — (A.3) 

sinhzlx — xq) 



or, after integration by parts 



S Q = —^ + 2(2 + ik) I dxe {2+ik)x — ( , 1 -. r. (A.4) 

u sinh 2x n / „ e 2(x-x ) _ e -2(x-* ) v ; 



Now, according to fl3.6|), if Co > °"i then xo > 0. Otherwise, it is necessary to adjust the 
background solution (|3.5|) by shifting xq through in/ 2 in order to be guaranteed that 
xq > 0. The singularity in the equation ( |3.5| ) is unimportant provided xq is positive. 
So, from now on it is assumed ctq > <j\. Therefore, xq is greater or equal to zero. But, 
x is less than zero, so < e^ x ~ x ^ < 1 and hence 



g2(x-x ) _ e -2(x-xo) 



_ e 2(x-x ) e in(x-x )^ (A. 5) 



n=0 

Substituting ( |A.5| ) in ( |A.4| ), we obtain 



i r0 00 

S = 2(2 + ik) e~ 2x ° / dx e^ +ik)x V e M*-*o)_ 

smh 2x J-oo ^ 

Clearly, the series is uniformly convergent so the above relation becomes 

oo „ 



(A.6) 



S = 2(2 + ik) e~ 2x ° V e" 4ra ° / rfx e (4+4n+ifc)x . (A. 7) 

smh2x ^ ./_«, 

After integration over x, we obtain 

So= — - +2i{2 + ik)e~ 2xo Y — ^- ■ — rr- (A.8) 

On the other hand, the above infinite series is a hypergeometric function. That is 

" e-4n, F(l,|fc+l,ifc + 2,e- 43; '») 

^ Q k-i(4 + 4n) k - Ai ^ ' ' 
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Therefore, we get the following relation 
So 



sinh 2xq 



The hypergeometric function is defined by |25| 



F(a,b,c,z) = f2 { ^^z n c + 0,-1, -2,... (A.ll) 



n=0 W 



where 

(a) n = v - , ; = a(a + l)...(o + w- 1) n= 1,2,3,... . (A.12) 
r(a) 

The above series defines a function which is analytic when \z\ < 1. Also, the derivatives 
of the hypergeometric function are given by 

d n (a) (b) 

—F(a, 6, c, z) = y Jn F(a + n,b + n, c + n, z). (A.13) 
dz n v 7 (c) n v 7 

Next, for n = 1 using ( |A.2| ) for the bulk three point coupling we have 

r° 

Si = 2 dx e (2+ ' k)x cosh 2{x - x ) coth 2{x - x ) (coth 2 2{x - x ) - l) . (A.14) 

On the other hand, if we differentiate the left hand side of ([A. 10 ) with respect to Xq, 
which is given by 



<99 r° / 

'" -/ dxe {2+ik)x Usmh2(x-x ) (coth 2 2(x - x ) - l) 

J — oo 

+8cosh2(x-x )coth2(x-x ) (l - coth 2 2(x - x ))) (A.15) 



dx 



and by comparing the above formula with ( |A.14j ), then the following equation may be 
derived 

Sl = + f° dxe^ 2+ik > —4 r- (A.16) 

4 ox J_ 00 smh2(x — x ) 

The second term in the above relation can be manipulated as before. Moreover, it is 
evident that we need to differentiate the right hand side of ( [A.10D which is equal to 

dSo 2 cosh 2x 

dx sinh 2 2x Q 

+4^ e ^F(l 1 iHl 1 ^ + 2 ) e"*») 
k — At 4 4 

+ 8 T^ e ~ 6X ° F ( 2 > 7 k + 2 ' 7 k + 3 ' e_4X °)- (A. 17) 

k — qi 4 4 
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Finally, by substituting the relation (A. 17) in ( |A.16| ), doing the computation of second 
term in the right-hand side of (|A.16|) and after simplifying we obtain 



cosh 2xo 

D\ = 2 

2sinh 2xq 



+ ^e- 2 -F(l,^ + l,^ + 2,e--o) 

+2^^ e~ 6x ° F(2, -k + 2 } -k + 3, e~ ix °). (A.18) 
k — 8i 4 4 



In the same way, we may derive ( A.l|) when n is equal to 2 or 3, however gradually 



the calculations become lengthy and we only write down the results, that is 

2 1 1 cosh 2 2x + 1 

02 — \ 1" q 

3 sinh2x 6 sinh 2x 

3 k-U 6 ni '4 fe + 1 '4 fe + ^ ] 

4f^| e ~" 0F(2 'I fc + 2 'I (A; + 3 ' e "" 0) 

-\ 7^t| F & \k + 3, Lk + 4, e" 4 -) (A.19) 



and 



13 cosh2x 1 cosh 2x + 5 cosh 2xo 
24 sinh 2 2x 24 sinh 4 2x 

4 1 f^ e ~ te ' F(2 '^ + 2 ^' + 3 ' e ~ 4 " ) 

+4-^^ e" 14 * F(4, ^ + 4, Ifc + 5, e~ 4a: °). (A.20) 
k — lot 4 4 

Appendix B 

In this Appendix we find the following integrals 

r° 

C n = rfxe jfcx cosh(v / 2/?0 o )coth n 2(x-x o ). (B.l) 

J — oo 

Here, cosh(v^2/30o) is proportional to the bulk four point coupling and is given by 

cosh(v / 2/50 o ) = (2 coth 2 2(x - x ) - l) . (B.2) 
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So, let us calculate such integrals 

o 

dxe ikx coth n 2(x-x ), (B.3) 

-oo 

where n — 1, 2, 6 . It is better to find the solution of the above integrals when n — 1. 
Considering the following inequality (see Appendix A) 

< e i{x - xo) < 1 

and therefore, in what follows we will use the expanded form of coth2(x — xq) as 

oo 

coth 2(x - x ) = 1 - 2 e in(x - x °\ (B.4) 

n=0 

It turns out to be simple if we consider this integral 

r° 

/ dxe (T+tk)x coth2(x-x ), (B.5) 

J — oo 

where r is a positive constant quantity which will be taken to zero at the end of the 
calculation. Moreover, by using ( |B.4| ), then ( |B.5| ) becomes 

/o t-0 °° 

dx e {r+ik)x _ 2 / dx e M*~*0) e (r+ik)x ( B £j 
■00 J-oc n=Q 

and regarding the arguments in Appendix A, we may evaluate the above relation to 
obtain 

00 

' + 2 y 1 — , % a ^ e ~ inx °- (b.7) 



k — ir k — (r + 4n)i 



Now, we can write down the desired result, that is, 

n=0 



On the other hand, the above series is equal to a hypergeometric function 



k- Ani k 4 4 

n=0 

and finally we find this formula 

j 1 = -i+|F(i,iife > iife+i )C -n (b.io) 

Now, let us compute ( |B.3|) when n = 2 and in order to solve it, it is sufficient to 
differentiate both sides of (|B.10| ) with respect to xq to obtain 

r_ - _L ^_ .-4-0 F ( 2 , -A; + 1, -k + 2, e- 4xo ). (B.ll) 



k k-Ai v '4 
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We can follow a similar method to get higher order forms of ( |B.3| ) which we need 
in this paper so, it is appropriate to write down all of them i.e. 

An <? <? 

' e~ 4xo F{2,-k+l,-k + 2,e- Axo ) 



k-Ai v '4 '4 

a-- j j 

-8x eVo _Uj_n o -4x ^ 



e -sx F ^ k + 2,-k + 3, e" 4 * ), (B.12) 



k-8i v 4 '4 



-4x 



16i 



F(2, -A; + l,-k + 2, e" 



-4x \ 



4 



k-8i 
16i 



8xo 



F(3, -A; + 2, -fc + 3, e' 
4 4 



-4xo^ 



(B.13) 



and 



o -4x p/o It II 1l i o -4x -> 



F(2,-A; + l, 7 A; + 2,e 



k-Ai y '4 '4 

^e~-oF(3,^ + 2,^ + 3,e--o) 



fc — 12i v 4 '4 

:>>2/ e- lfo0 F(5,-fc + 4,-fc + 5,e- 41 ' ) (B.14) 



16i v 4 '4 



1 2? 7 7 

e" 4:co F(2,-fc + l,-A; + 2,e- 4xo ) 



fc A;-4z v '4 '4 

e" 8:E0 F(3, -A; + 2, -fc + 3, e" 4 ^; 



fc-8i v 4 '4 
1121 e~ 12xo F(A -k + 3 -k + A e~ ix °) 

-^e--F(5,^ + 4,l, + 5,e--) 

- , 64 ' e- 2Cte0 F(6, -A; + 5, -jfe + 6, e" 4 * ). (B.15) 
— 20i K A '4 ' ; v ' 
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